Characterizing Charge Diffusion in CCDs with X-rays 
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S : ABSTRACT 

■^ . We demonstrate the effectiveness of two techniques for using x-rays to evaluate 

p ! the amount of charge diffusion in charge coupled devices (CCDs). We quantify 

O 



the degree of charge diffusion with two parameters: (Jd-, the standard deviation 



5-1 I for a Gaussian diffusion model, and Q, a ratio of the point spread function (PSF) 

^ ■ peak to its wings, ad and Q are determined by fitting a model to a pixel energy 

histogram, and by summing the PSF of all x-ray events, respectively. Using seven 
test devices, we investigate the precision of these two techniques and demonstrate 
rS \ that they produce compatible results. The histogram fitting method is sensitive 

d ' to the structure of the electric field within these devices, in addition to the 

inherent charge diffusion properties. The Q ratio is a very simple parameter to 
measure and provides an easily accessible method for quickly evaluating a CCD's 
diffusion length. 

Subject headings: : instrumentation: detectors - methods: laboratory - diffusion 

1. Introduction devices it is extremely important to have 

a set of standard methods for testing and 

The use of Charge Coupled Devices quantifying their properties. There are a 

(CCDs) as light collectors for photometry number of useful CCD characteristics that 

and spectroscopy is ubiquitous in modern ^^^ relatively straightforward to measure, 

observational astronomy. To understand p^^ example, carefully refined techniques 

fully the limits and capabilities of these fo^ ^gg^i^g ^^e gain, linearity, and charge 
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Fig. 1. — Cartoon representation of an x-ray 
absorption in a CCD and the resulting charge 
diffusion cloud. An x-ray enters through the 
backside surface at the top, and interacts at 
some height z above the gates. The initial 
size of the electron cloud is small compared 
to the size of the pixels. As the electrons are 
pulled downward at a velocity proportional to 
the electric field strength they random walk 
horizontally. A cross-section of the cloud at 
any height z in the silicon looks approximately 
like a two dimensional Gaussian distribution 
of electrons (insets 1, 2, and 3). When the 
cloud reaches the bottom of the device it is 
segmented by the grid of the pixel structure. 

transfer efficiency (CTE) have been de- 
veloped (Janesick 2001). One important 
CCD trait for which an efficient and prac- 
tical testing procedure does not exist is 
charge diffusion. 

Figure 1 depicts schematically the 
process of charge creation and collection 
within a CCD. When a photon reaches the 
surface of a CCD it penetrates through 
some depth of silicon before it is absorbed. 
At the point of interaction this photon lib- 
erates one or more charge carriers. A pho- 



ton with energy less than ~4 eV creates 
one charge carrier. Higher energy x-rays 
produce a packet of charge in which the 
number of charge carriers is proportional 
to the energy of the photon. These charge 
carriers are either electrons or holes, de- 
pending on the CCD design. Hereafter we 
will assume that the buried channel is in an 
n-type semiconductor, so the charge car- 
riers are electrons. The electron thermal 
velocities cause them to wander away from 
their original location into neighboring pix- 
els. This spreading effect is what is termed 
charge diffusion. The electrons produced 
by photon absorption are eventually con- 
verted into a digital image, so charge diffu- 
sion places a fundamental physical limit on 
the resolution in the final image produced 
by a CCD. 

Kraft et al. (1995) characterized charge 
diffusion in the process of evaluating quan- 
tum efficiency in x-ray CCDs, but their 
method involved a synchrotron source to 
produce a continuum of x-ray energies. 
Janesick (2001) also describes several tech- 
niques, such as measuring the modulation 
transfer function by projecting a sinusoidal 
light pattern across the CCD. In general, 
these methods can be effective, but they 
require either expensive specialized equip- 
ment or a significant investment of time. 
Such procedures are better suited to CCD 
manufacturers or engineering labs. De- 
spite the demonstrated importance of un- 
derstanding charge diffusion effects when 
designing astronomical instruments, there 
is currently no method for effectively eval- 
uating charge diffusion that is accessible to 
the astronomical community without spe- 
cialized equipment. 

The techniques presented here rely on 



the use of soft x-rays from a mild radioac- 
tive source as a means of producing dis- 
crete, single photon events. These events 
are approximately point sources, since the 
size of the initial charge cloud is generally 
very small compared to the size of the pix- 
els (Tsunemi et al. 1999). Of course, the 
experimenter has no control over where the 
x-rays land on the CCD surface or at what 
depth in the silicon each photon produces 
charge, so the x-ray events are randomly 
distributed throughout the volume of the 
CCD. Thus, we can gain little informa- 
tion by fitting a model to an individual 
event. An x-ray interacting near the sur- 
face of a low-diffusion device will spread 
charge over a large area, much like an x- 
ray interacting deeper down but in a de- 
vice with higher charge diffusion. We can, 
however, break this degeneracy by model- 
ing the entire ensemble of x-ray events and 
including the absorption probability as a 
function of depth (i.e. absorption proba- 
bility increases exponentially with depth, 
so there will be far more interactions deep 
in the silicon than near the surface). For 
the two methods described here, we com- 
pare the histogram of pixel energies or the 
stacked and summed x-ray PSF to models 
that incorporate all of the available infor- 
mation about the charge collection process. 
In the following section we describe the 
experimental setup and data collection 
procedures, as well as the formula for the 
calculation of the Q ratio. §3 details the 
creation and interpretation of K-a pixel 
histograms, while §4 explains the theoreti- 
cal model used to fit those histograms. We 
present results for the two methods in §5 
and final conclusions in 56. 



2. Experimental Setup and Proce- 
dures 

Seven devices, labeled A-G, were used 
for the diffusion measurements. The first 
three are CCID-20 chips, chosen primar- 
ily because of their disposability - these 
CCDs were not likely to be used in a 
science-grade instrument, so damages in- 
curred while testing them would not have 
been disastrous. The other four are CCID- 
45 Orthogonal Transfer Array (OTA) de- 
vices (Tonry et al. 1997). Each OTA has 
64 individually addressable cells of approx- 
imately 500x500 pixels each. Table 1 lists 
the relevant information for all of the de- 
vices used. This set of seven CCDs pro- 
vides a useful variety of physical param- 
eters that affect the measured charge dif- 
fusion, such as silicon thickness and pixel 
size. The OTAs have the unique additional 
feature of allowing the voltage level of the 
substrate to be adjusted during testing, 
which provides another way to modify the 
diffusion properties while holding all other 
parameters constant. 

To provide soft x-ray illumination, we 
use a 1 mCi Fe55 x-ray source. The Fe55 
isotope emits x-rays when the nucleus cap- 
tures an inner electron and becomes man- 
ganese. The most important x-ray emis- 
sions are the Mn K-a event at 5.899 keV 
and the K-/9 line at 6.490 keV. K-a and 
K-/9 x-rays interact in the silicon and pro- 
duce a characteristic charge cloud equiv- 
alent to 1620 and 1778 electrons, respec- 
tively. Each of the four devices is mounted 
in a dewar that has a beryllium window 
in its faceplate, opaque to optical and 
infrared light but transparent to x-rays. 
By placing the Fe55 x-ray source directly 
against this window while the CCD is inte- 



Table 1: Devices Tested 



Device Model Thick. Pix. Size 
(CCD-) (CCID-) (/xm) (fim) 



hbf* Resistivity Notes 
(pm) (kn — cm) 



A 


20 


20 


15 


4 


0.3 


B 


20 


45 


15 


6 


5 


C 


20 


45 


15 


8 


5 


D 


45 


45 


10 


10 


5 


E 


45 


45 


12 


10 


5 


F 


45 


45 


12 


12 


14 


G 


45 


75 


12 


14 


14 



bright defects; useable area ~oo7o 
bad parallel CTE along one side; 
useable area ~90% 
tree-ring pattern fab. error 
variable substrate voltage 
variable substrate voltage 
variable substrate voltage 
variable substrate voltage 



'From histogram fits. See section 4 



grating, we get a dark image with hundreds 
of isolated single x-ray events. 

Figure 2 shows a simple histogram of all 
pixel values from an x-ray image of CCD- 
A. This histogram includes every pixel 
in the image (above a threshold of 1000 
ADU), regardless of whether it is identified 
as a K-a pixel or not. Such histograms are 
useful for accurately measuring the gain 
of a device, but to measure charge diffu- 
sion we really want to limit ourselves to 
those pixels in the immediate vicinity of 
an x-ray event. To produce the required 
x-ray pixel histogram, an x-ray image is 
scanned with a point source locating al- 
gorithm, which locates all pixels that are 
adjacent to a K-a x-ray event. This algo- 
rithm first creates a raw histogram, such 
as the one in Figure 2, and locates the K- 
a peak to determine N^a, the number of 
ADUs that corresponds to a complete K-a 
charge cloud (i.e. 1620 electrons times the 
gain of the device). We then examine all 



pixels that have a signal more than seven 
times the RMS noise of the background, 
along with their eight immediately adja- 
cent pixels. For each patch of nine pixels, 
we subtract a constant background level 
and calculate the sum of nine pixel values, 
N patch- If this sum is close to A^fco, meaning 
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then all nine pixels are tagged as 'K-a pix- 
els' and are included in the calculation of 
Q and in the x-ray pixel histogram. 

To measure Q, the 3x3 arrays of K-a 
pixels, centered on the brightest pixel, are 
stacked and summed to produce a compos- 
ite x-ray PSF. Q is then defined as the nor- 
malized ratio of the sum of the eight outer 
pixel values to the peak pixel value in the 
composite PSF: 



Q 



1 ^ {eight outer pixels) 
8 {central pixel) 




Fig. 2. — The complete pixel value histogram 
for a typical x-ray image. The Mn K-a and 
K-/3 peaks are known to produce 1620 and 
1778 electrons, respectively, so the gain of this 
device can be easily determined by measuring 
the peak location in ADU and dividing the 
appropriate electron number by this value. 

Thus, for a device with zero charge diffu- 
sion — so that all 1620 electrons produced 
by a K-a x-ray are collected in a single 
pixel every time — the Q ratio is zero. For 
infinite charge diffusion — so that each of 
the nine pixels receives an equal fraction of 
the total charge — Q is unity. 

3. The X-ray Pixel Histograms 

After computing Q, the algorithm then 
collects the complete set of all K-a pixels 
(meaning all nine pixels from every patch 
centered on a K-a event) and bins them 
into a histogram indicating what fraction 
of the x-ray's energy was deposited in each 
pixel. To facilitate comparison with other 
devices and models, these x-ray pixel his- 
tograms are normalized such that the inte- 
gral over all histogram bins is unity. These 
histograms are dominated by the back- 
ground pixels, clustered around zero en- 
ergy, so we display them on a log scale 



-2 


-1.0 


i«; 




Qh 


-1.5 


5^ 


-7 












-2.3 


3 
^ 


-3.0 


o 


-3.5 


O 






-4.0 


"^ 


-3 










a. 


-3.2 


"^ 




o 


-3 4 


>^ 




^ 


-3 6 


s 




fe 


-3.8 


o 


-4.0 


o 






-4.2 





— CCD-D:45/10 

— CCD-F: 45 / 12 

— CCD-G: 75/12 . 




t^- 




0.6 0.7 0.8 0.9 1.0 
Fraction of k-alpha energy 



Fig. 3. — Logarithmic x-ray histograms from 
CCD-D, F and G in red, black and blue, 
respectively. All three histograms were cre- 
ated with the CCDs set to the same sub- 
strate voltage of -10 V. The legend in the up- 
per panel gives the thickness and pixel size 
in nm for each device. CCD-D and F have 
the same thickness, but different pixel sizes. 
CCD-D, with smaller pixels, shows a dearth of 
events with highly concentrated charge collec- 
tion. CCD-F and G have the same pixel size, 
but different thicknesses. The thicker device, 
CCD-G, has substantially fewer x-ray events 
with more than 70% of the energy collected by 
a single pixel, relative to the thinner CCD-F. 



to highlight the subtle features near unity, 
which are strongly affected by the amount 
of charge diffusion. Figure 3 shows three x- 
ray histograms, demonstrating how phys- 
ical properties of the CCD alter the his- 
togram appearance. 

The prominent features of these x-ray 



histograms can be explained in terms of 
the physical process of charge diffusion and 
collection. The strong peak around an 
energy value of zero is of course due to 
the pixels included in each 3x3 patch that 
have very little or no contribution from 
the nearby K-a electrons. This peak ex- 
tends into negative energy values because 
of noise fluctuations that leave some pixels 
below the median background level. To the 
right of this peak, as should be expected, 
the slope of the histogram is mildly nega- 
tive through the intermediate energy bins. 
The first feature of interest to the deter- 
mination of charge diffusion is the slope 
discontinuity, which is particularly promi- 
nent in CCD-G as shown in Figure 3.^ This 
slope discontinuity indicates a sudden de- 
crease in the number of pixels containing 
a large fraction of the K-a electrons. The 
energy at which this critical point occurs 
marks the energy remaining in the central 
pixel when an x-ray lands at the center of 
a pixel and the surface of a CCD. In any 
CCD, there are essentially two ways for an 
x-ray event to deposit a large fraction of 
its charge in a single pixel: it can land very 
near the center of the pixel's square cross- 
section (allowing it lots of room to spread 
out within a single pixel), or it can interact 
with the silicon at a vertical height close 
to the gates (giving it very little time to 
spread out to other pixels). In CCDs with 
small pixels (such as CCD-D in Fig. 3) the 
first avenue is less effective, so the disconti- 
nuity point is translated to the left relative 
to large-pixel devices. In thick devices, the 
second effect is suppressed because the x- 



^We emphasize here that these histograms are plot- 
ted on a logarithmic scale, so these are in fact very 
subtle slope discontinuities. 



rays must travel through a greater volume 
of silicon to reach the requisite depths for 
containing the charge cloud in one pixel. 
Thus, thicker devices (such as CCD-G in 
Fig. 3), have fewer single pixels holding 
most of the K-a electrons. 

The final peak of these x-ray histograms 
is populated by single-pixel x-ray events. 
These are events in which the entire charge 
cloud is collected in one pixel. The pres- 
ence of this peak in all of our x-ray his- 
tograms indicates that there is a sizeable 
probability for an x-ray penetrating deep 
into the CCD to deposit all of its energy 
into one pixel, rather than to spread a 
small fraction of it to neighboring pixels. 
As explained in §4, this trait is amplified by 
the presence of pixel barrier fields, which 
effectively focus charge clouds into pixel 
wells. 

4. The Diffusion Model 

To construct a model for comparison 
with the histograms described above, we 
consider an x-ray incident on the backside 
of a CCD that is absorbed at some height 
z above the frontside gate, as shown in 
Figure 1. The potential difference across 
the device from the backside to the buried 
channel produces an electric field that 
draws the electrons into the potential well 
which defines each pixel. If these were free 
electrons in empty space, they would be 
accelerated by this electric field. Inside a 
silicon lattice, however, these charges are 
impeded by collisions with lattice atoms 
and it can be shown that the drift velocity 
is proportional to the electric field at that 
height: v^ oc E{z). (Janesick 2001, p.478) 
As the electrons are gradually pulled down- 
ward, they also execute a random walk 



horizontally, so that when they are finally 
collected in the gate structure the result 
is a substantially more extended profile 
than the original charge cloud. This final 
charge distribution is generally taken to 
be a Gaussian (Hopkinson 1987; Janesick 
2001), although some calculations have 
shown that it should be more centrally 
peaked (Pavlov & Nousek 1999; McCarthy 
et al. 1995). For simplicity, we restrict 
ourselves to the classic Gaussian charge 
distribution in all models. 

We can take the RMS horizontal dis- 
placement to be proportional to the square 
root of the time it takes to reach the 
frontside gate at z=0: 
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In the simplest case, we can take the elec- 
tric field to be constant with depth, so that 
a device with thickness d and potential dif- 
ference V from frontside to backside has 

E{z) = V/d and we can now write 
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where a^ is the standard deviation of the 
Gaussian diffusion cloud that would result 
from an x-ray interacting at the surface, 
where z = d. Now for a given value of a^ 
and a K-a x-ray event at a known height z 
we can create a final Gaussian distribution 
of charge. If we also choose an (x,y) posi- 
tion for the center of the cloud relative to 
the pixel boundary, we can integrate that 
final charge distribution over a 3x3 array 
of CCD pixel boundaries. To create a his- 
togram from these pixel values, we want 



to repeat this process for a large number 
of x-ray events spread over an appropriate 
range of (x,y,z) values. For this, we need 
to know what fraction of x-ray events is 
likely to occur at any given depth in the 
silicon. K-a photons in silicon have an e- 
folding absorption depth of 27.8 /xm, so in 
a device of thickness d /i?T2 the fraction of 
absorption events occuring between height 
Z\ and Z2 from the bottom is: 
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(4) 



The last physical effect that we must 
consider is the presence of barrier fields 
produced by the gate structure on the 
frontside of each device. The division of 
any CCD into pixels necessarily requires 
potential barriers which delineate the in- 
dividual pixel wells. These barriers intro- 
duce a transverse electric field that pre- 
vents electrons very deep in the silicon 
from crossing over into neighboring pix- 
els. We simulate this effect by introduc- 
ing a barrier field height parameter, /i^/, 
such that once an electron reaches a height 
z < hbf it can no longer diffuse away from 
its current position. This effectively trun- 
cates the diffusion process at a height h^f 
above the buried channel. To accommo- 
date this, we rewrite Eq. 3 as: 
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(5) 

Thus, in our models, a 75 /xm thick device 
with a barrier field height of 11 fim exhibits 
similar diffusion characteristics to a 65 fim 
thick device with a 1 fim barrier field. The 
one important difference is that the final 11 
fim, of the first device will absorb substan- 



tially more x-rays than the final 1 ^m of 
the second device, meaning that there will 
be more single-pixel x-ray events populat- 
ing the final peak of the x-ray histogram. 
Hence, the effect of increasing /i^j in our 
models is to amplify the final histogram 
peak at the expense of the energy bins just 
below the peak. 

To produce a model energy histogram, 
we divide the silicon thickness into thin 
slabs. For each slab, Eq. 5 determines 
the apparent (t(z) at that depth, and Eq. 
4 sets the relative number of events that 
have that cr(z). Next, we create a regu- 
lar grid of (x,y) values arranged across one 
triangular eighth of the surface of a single 
pixel. Due to the symmetry of a square 
pixel, this is effectively the same as us- 
ing a grid that extends across the entire 
pixel surface. For each (x,y) pair, we use 
the a{z) appropriate to the current slab 
to create a final Gaussian diffusion cloud 
projected onto the pixel grid. We then in- 
tegrate over each pixel to determine the 
total x-ray counts detected for each event, 
and collect these pixel values into a nor- 
malized histogram in exactly the way that 
was done for the actual data. The final 
result is an x-ray histogram model with a 
known ad that can be directly compared to 
the real histogram. 

To fit a model to each device we re- 
peat this process for a wide range of aa 
and hhf. The silicon thickness and intrin- 
sic noise are fixed parameters for each de- 
vice. From the resulting library of model 
histograms, we choose the curve that most 
closely matches the data histogram and 
record that model's value of ad as our mea- 
sured a hist- Figure 4 shows an actual his- 
togram for CCD-G along with the three 




Fig. 4. — Demonstration of the precision of 
our model histogram fitting technique. An ac- 
tual x-ray pixel value histogram from CCD-G, 
operated with a substrate potential set to - 
20 V, is shown in the solid step-line. Three 
model histograms are overplotted as dotted, 
dashed, and dash-dot lines. These models 
have ad values of 0.28, 0.30, and 0.32, respec- 
tively, in units of the pixel width. For this 12 
fim pixel device, these correspond to 3.36, 3.6, 
and 3.84 fim. All three models have a barrier 
field height of /i^j = lA^m 

closest matches from our model histogram 
library. The range in ad for these three 
models is 4% of the pixel width, or 0.48 /xm 
for this 12 /im device. From these model 
histogram fits we can determine that the 
precision of a^ist is ±0.24/im. All three 
models use a barrier field height of 14 ^m, 
which is 19% of the total 75 /im thickness, 
and is close to the pixel width. In general, 
all of the best fitting models for a given de- 
vice agree to within 1 ^m on the height of 
the barrier field, regardless of the substrate 
voltage. 

The measured barrier field heights for 
these seven devices are of the order of the 
size of the pixels. As a general rule, chang- 
ing the barrier field height in our model 
has no effect on the position of the slope 



discontinuity point, which is the primary 
indicator of ad. The height of the final 
peak and the depth of the preceding trough 
in each histogram are the sole indicators 
of hbf. In devices with very low charge 
diffusion, such as the 20 /im thick CCD- 
A, our determination of h^f is less pre- 
cise, because x-ray events have to be very 
close to the pixel boundaries for the bar- 
rier field to have any effect due to the in- 
herently small size of the diffusion cloud. 
These low-diffusion devices will also yield 
somewhat less accurate measures of aust- 
When there is very little charge diffusion 
the slope discontinuity point that indicates 
the value of ahist is very close to the final 
peak in the histogram, and can therefore 
be confused by the presence of the preced- 
ing trough feature. Table 1 lists the barrier 
field heights measured from histogram fit- 
ting in column 5. 

For the Q measurements, we simply 
need to translate our measured Q ratio into 
a physical measurement that can be di- 
rectly compared to the (Thist value in mi- 
crons. For this purpose we employ the 
same setup to produce a model aggregate 
PSF, and perform the sum and division 
of the nine pixels to get a value of Q. 
For each device the three physical prop- 
erties of thickness, pixel size, and barrier 
field height are used as parameters for the 
model to define a curve that shows the in- 
put value of ad plotted against the output 
Q. Figure 5 shows these curves for three 
of the seven devices tested. From this plot 
we can take our unitless measured Q value 
and translate it directly into a value of aq 
in ^m. 
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Fig. 5. — The results of our model recreation 
of the Q calculation process for three of the 
seven CCDs. For each CCD, the thickness, 
pixel size and barrier field height are set as 
fixed parameters. Then a series of Ud values 
are input, and our model creates a summed 
PSF made up of nine pixels. This model 
summed PSF is evaluated to get Q, which is 
then plotted here against the input value of 
Gd- These curves are used to translate our 
measured values of Q into a corresponding aq 
in ^m. 

5. Results 

Using the procedures described above, 
we determine a value of ahist from his- 
togram fitting, and a value of aq from the 
summed PSF for each device (and at mul- 
tiple substrate voltages in the case of the 
OTAs). These results are shown in Fig- 
ure 6. The upper panel of Figure 6 plots 
the ahist values determined from histogram 
fitting against the substrate voltage for all 
seven devices. The three non-OTA devices 
(CCD-A, B and C), are plotted as single 
points at Vsuh = 0. The lower panel of Fig- 
ure 6 provides the same plots for aq as a 
function of substrate voltage. 

From these plots we can discern a num- 
ber of important effects. First, our two 
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Fig. 7. — As in Figure 6, but showing only gq 
for CCD-D and CCD-G. The model of Hop- 
kinson (1987) (see Equations 6 and 7) is over- 
plotted for each set as a solid line. For both 
devices we set Vgate = 8.51/ and we use a dop- 
ing concentration oINa = 8.9 x 10^^ for CCD- 
G and Na = 2.5 x lO^^ foj. cCD-D, based on 
the resistivities shown in Table 1. 



Fig. 6. — (TQ and atist plotted as a function 
of substrate voltage. The non-OTA CCID-20 
devices are plotted as single points along Vsub 
= V. A single letter on the right side of the 
figure labels each curve with the correspond- 
ing device name. Typical error bars are shown 
in the upper left corners. 

methods both verify that thickness has a 
strong affect on a device's charge diffusion. 
The 75 nm thick CCD-G has almost a fac- 
tor of two more charge diffusion than the 
45 /xm OTA devices D, E, and F. CCD- 
A, at a thickness of 20 fim has more than 
a factor of two less diffusion than the 45 
fj,m CCD-B and C. Secondly, the distinct 
histogram shapes resulting from differences 
in pixel size that are shown in Figure 3 
produce the same physical diffusion length 
when we take into account the pixel size 
in our translation to 0"^ in fim. Thus, the 
curves for CCD-D and F are almost identi- 



cal to within the error bars. From the four 
OTAs we can evaluate the effects of chang- 
ing the substrate voltage. As expected, 
a more negative substrate bias establishes 
a stronger electric field, which results in 
lower measured values of the charge diffu- 
sion. This is true for all devices across the 
whole range of sampled voltages (except 
for one erratic point on CCD-E). 

Hopkinson (1987) provides a theoreti- 
cal model for the size of a charge diffu- 
sion cloud, 0", as a function of the operat- 
ing voltage, V. This model starts from the 
assumption that the electron drift velocity 
is linearly proportional to the electric field 
strength at any height in the silicon — an 
assumption we also adopt in Equation 2. 
Hopkinson (1987) predicts the final la ra- 
dius of a charge cloud that begins at the 
surface of a CCD as: 
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where k is Boltzmann's constant, T is the 
temperature of the CCD, e is the elec- 
tric permittivity of sihcon (1.044 x 10~^^ F 
cm~^), e is the electron charge, and d is the 
thickness of silicon. A^^ is the concentra- 
tion of acceptor atoms, which we determine 
from the resistivity values quoted in Table 
1 using equation 1.24 from Janesick (2001, 
p. 75). dd is the depletion depth, given by: 
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where V is the operating voltage of the de- 
vice (approximately the gate voltage mi- 
nus the substrate voltage). Figure 7 shows 
the data for CCD-D and CCD-G with 
the above model overplotted as solid lines. 
This figure demonstrates that the observed 
variation with Vsuh and d is reasonably well 
matched by the Hopkinson model for large 
values of the substrate voltage. When Vsub 
is small, the devices approach underde- 
pletion, so that the depletion depth, dd, 
is very close to the device thickness, d. 
Hopkinson's model does not apply in this 
regime, which is why the data begin to de- 
viate from the theory. 

When comparing the two methods 
against each other, we find that they pro- 
duce consistent, if not identical results. 
Figure 8 plots the measured value of ahist 
against aq for all seven devices. If the 
two methods measured precisely the same 
value for the physical diffusion length, the 
points would fall along the solid line in this 
figure with a slope of one. A least squares 
fit to the points reveals a linear relation of 
the form ahist = 0.83 x aq + 0.82, shown 
as the dashed line. As noted above, we 
have good reason to be suspicious of the 
accuracy of the histogram fitting model 
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Fig. 8. — aq plotted against ahist for all seven 
devices. The solid line shows where the points 
would fall if the two methods were in perfect 
agreement. A least squares fit is plotted as the 
dashed line, with the functional form given in 
the lower right. 

at very low values of ad- When ad is 
less than about 0.25, the prominent his- 
togram features become crowded, leading 
to an inaccurate estimate of the diffusion 
length. This deficiency could be corrected 
by employing substantially more x-rays in 
the histogram to improve energy resolution 
and signal to noise. 

6. Conclusions 

Using a set of four devices with a range 
of thicknesses and charge collection char- 
acteristics, we have evaluated two meth- 
ods for using x-rays to measure the amount 
of charge diffusion in CCDs. With a rela- 
tively simple model for the charge diffusion 
effect built off of the assumption of a Gaus- 
sian diffusion cloud, we are able to produce 
model histograms that match the observed 
pixel energy distributions very accurately. 
By varying the input parameters of the 
model to fit the observed data we can put 
precise constraints on both the amount of 
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charge diffusion and the height of the bar- 
rier fields around each pixel. We introduce 
the Q ratio to measure the amount of dif- 
fusion based on a simple quotient of the 
outer eight pixels to the central pixel in 
the summed x-ray PSF. This method pro- 
vides an even more direct means of assess- 
ing the diffusion properties of each device, 
as it requires no fitting procedure and is 
independent of model assumptions for rel- 
ative diffusion measurements. Using the 
same model as for the histogram fitting 
method, we are able to turn this simple 
Q ratio into an absolute measurement of 
the size of the diffusion cloud. Comparing 
the fundamental diffusion parameters from 
each method, cThist and ctq, shows that the 
two methods produce results that are in 
good agreement. Due to the confusion of 
histogram features in low diffusion devices, 
ahist is best for measuring relatively large 
values of a^ and ctq is more appropriate 
for the smaller values. X-ray images are 
easy to collect in the lab, requiring very 
little specialized equipment, and only min- 
imal image processing is necessary to get 
the data needed for these techniques. All 
of these advantages make the x-ray his- 
togram and PSF ratio methods very useful 
"quick-look" additions to the CCD evalu- 
ation toolkit. 
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